WHOI-93-46 


Bottom  Boundary  Layer  Stress  Measurements 
with  BASS  Tripods: 

Data  Report  STRESS  1988-89 

nett  » 

Thomas  F.  Gross  and  Albert  J.  Williams  3rd 


A 


Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Massachusetts  02643 


November  1993 

Technical  Report 


Funding  was  provided  by  the  Office  of  Naval  Research  under  Contract 
Nos.  NO0014-89-J-1068  and  N00014-90-J-1046. 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United  States 
Government.  This  report  should  be  dted  as  Woods  Hole  Oceanog.  Inst.  Tech.  Rapt., 

WHOI-93-46. 


Approved  for  public  release;  distribution  unlimited. 


Approved  for  Distribution} 


Accesion  For  ) 

NTIS  CRA&I 

OTIC  TAB 

□ 

Unannounced 

□ 

Justification 

J 

Authors*  Note:  In  this  report,  die  height  above  bottom  of  the  top-most  sensor 
of  BayShoct  IV  is  specified  as  440  cm  above  bed.  Subsequent  to  the  generation 
of  this  report,  it  has  been  noted  that  the  field  log  indicated  the  height  above  bot¬ 
tom  as  499  cm  above  bed  The  499  cm  height  has  been  verified  and  should  be 
corrected  by  the  reader  on  Page  2  and  in  Figure  1. 

July  1993 


ii 


page  1 


Bottom  Boundary  Layer  Stress  Measurements  with  BASS  Tripods 
Data  Report  STRESS  1988-89 

Thomas  F.  Grossf 
A.  J.  Williams  Olft 


Overview 

The  Benthic  Acoustic  Stress  Sensor  tripods,  BASS,  are  bottom  landing 
instrument  platforms  designed  to  measure  the  turbulent  bottom  boundary  layers. 
The  turbulent  bottom  boundary  layer  is  the  turbulent  flow  in  the  region  from 
the  sediment  bed  to  5  to  10  meters  above  die  bottom.  Within  this  region  tur¬ 
bulent  stress  slows  the  flow  from  midwater  column  speeds  to  zero  at  the  "no- 
slip"  condition  of  the  sea  bed.  Boundary  shear  stress  scales  sediment  erosion 
rates  and  provides  boundary  friction  which  affects  the  full  water  column 
momentum  balance.  The  stress  is  due  to  the  mean  boundary  layer  shear  and 
enhanced  by  wave  generated  turbulence.  Therefore  measures  of  the  mean 
profile  and  wave  action  are  needed,  as  well  as  a  modeling  method  by  which  the 
bed  shear  stress  can  be  derived.  The  Grant,  Madsen  and  Glenn  models  of 
wave-current  interaction  will  be  used  for  this. 

The  sheared  velocity  profile  is  to  first  order  a  logarithmic  profile  described 
by  U(zy>*ufKlnz/2c.  The  shear  velocity  scale,  iu,  describes  the  turbulence 
intensity  and  is  used  to  derive  die  boundary  shear  stress  x »=pu»2.  These  quan¬ 
tities  are  measured  by  obtaining  velocity  profiles  at  half  hourly  intervals 
through  the  bottom  5  meters  of  the  flow. 

The  wave  action  is  derived  by  obtaining  pressure  power  spectra  every  half 
hour.  Bom  this  a  peak  frequency  and  pressure  variance  are  obtained.  By 
using  linear  wave  theory  these  are  related  to  the  near  bed  wave  velocity  and 
excursion  amplitude.  The  validity  of  these  conversions  are  checked  by  com¬ 
parison  to  velocity  variance  which  is  dominated  by  wave  motions. 

The  BASS  tripods  provide  turbulent  bottom  boundary  layer  descriptions 
by  measuring  the  velocity  at  six  fixed  locations  within  5  meters  of  the  bed. 
The  mean  velocity  data  is  used  to  derive  logarithmic  profile  estimates  of  n».  In 
addition  die  BASS  velocity  data,  which  is  sampled  at  1.6  Hz,  can  provide  esti¬ 
mates  of  kinetic  energy  and  die  Reynolds  stress  tensor  by  averaging  velocity 
component  cross  correlations: 

1  Xmin  _  _ 

<“>>-3—  I  UiWjiOdt  -  U,Uj 

where  Uj  (1 )  j,  die  original  time  series  data.  Ui  is  the  avenge  over  30  minutes 
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and  u'j  is  the  turbulent  velocity  component. 

Data  Results 

The  BASS  tripods  were  deployed  twice  at  the  STRESS  site.  Table  1. 
During  the  first  deployment  one  six  BASS  tripod  was  put  down  at  the  90m  C3 
station,  data  tape  T2000.  The  second  deployment  put  that  same  tripod  at  the 
90m  site,  data  tape  T2006,  and  another  four  BASS  tripod  at  the  inner  shelf  C2 
station  in  SSm,  data  tape  T2005.  Calibration  zeros  were  obtained  for  both  tri¬ 
pods  after  the  second  deployment  Zeros  were  obtained  by  bagging  the  instru¬ 
ments  to  prevent  flow  through  the  sensors  and  placing  them  on  the  bottom  at 
the  Sausalito  dock. 

BASS  tripods  were  deployed  at  the  STRESS  sites  winter  1988-89.  Three 
data  tapes  were  recovered: 

T2000: 

Turn  on:  Nov.  28  18:30  GMT  =  332.7708  year  day  (y.d.) 

In  water.  Nov.  28  21:30  GMT  1988  *  332.8958  y.d. 

to  Jan.  24  9:30  GMT  1989 
90  meter  site,  Bayshore  IV  BASS,  Clean  tripod. 

Six  BASS  pods  at  22.6,  53.1, 113.2,  198.9,  259.4,  440  cm  above  bed. 

Two  Transmissometers  at  46  and  144  cm  above  bed 
Eight  thermistors  at  40,  65, 113, 204, .... 

Pressure  sensor  4  cm  above  bed,  serial  number  13032 
Compass:  unknown,  not  recorded 

T2006: 

Turn  on:  Jan.  26  21:00  GMT  1989  »  25.8750  y.d.1989  =  390.8750  y.d.1988 
to  Mar.  18  15:30  GMT,  1989 
90  meter  site,  Bayshore  IV  BASS,  Clean  tripod. 

Six  BASS  pods  at  22.6,  53.1, 113.2,  198.9,  259.4, 440  cm  above  bed. 

Two  Transmissometers  at  46  and  144  cm  above  bed 
Eight  thermistors  at  40, 65, 113,  204, .... 

Pressure  sensor  4  cm  above  bed,  serial  number  13032 

Compass:  266°,  Rotate  U  to  EAST  with  addition  of  90-266-25-16.9=152° 

T2005: 

Turn  on:  Jan.  25  19:00  GMT  1989  =  24.7917  y.d.1989  =  389.7917  y.d.  1988 
In  water.  Jan.  25  22:00  GMT  1989  *  24.9167  y.d.1989 
to  Mar.  17  22:30  GMT  1989 
55  meter  site.  Glass  BASS,  LDV  tripod. 

Four  BASS  pods  at  41.5,  76.5, 111.5,  510.  cm  above  bed. 

Two  Transmissometers  at  50  and  148  cm  above  bed 
Three  thermistors  at  46,  112 , 510  cm  above  bed 
Pressure  sensor  4  cm  above  bed,  serial  number  9126 
Compass:  unknown,  compass  defective 
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Data  recorded  on  each  tape: 

Month  day  hour  minute 

Half  hour  averages  of  1.6  Hz  velocities  along  each  of  four  axes  for  each  pod. 
Quality  word  for  each  axes  of  each  pod. 

Half  hour  average  cross  correlations  of  1.6  Hz  data  of  pod  axes:  aV,  a’b\  a’c\..etc. 
Ten  minute  readings  of  thermistors  and  transmissometers 
Two  minute  averages  of  U  and  V  at  one  level 

Fourier  transformations  of  pressure  variance,  frequencies  (l:40y307.2s  .'fairly  good. 
Fourier  transformations  of  velocity  component  variance, 
frequencies  (l:40)/307.2s  problematic. 

On  the  second  deployment  only  (T2006,  T2005)  the  following  were  added: 

Compass 

Half  hourly  pressure.  Not  averaged,  so  terribly  aliased  by  surface  waves. 


Calibrations: 

BASS  velocities 

BASS  velocity  calibrations  involve  a  gain  factor  and  a  zero  offset  The 
gain  factor  is  only  a  function  of  the  path  length  and  speed  of  sound  in  water, 
neither  of  which  are  variable,  gain  *  0.00368  cm/sec  /bit  The  zero  offsets  are 
due  to  capacitence  variations  of  the  cables.  This  changes  slightly  with  each 
configuration  of  the  tripod.  Therefore  aero  velocities  are  measured  by  deploy¬ 
ing  the  tripod  with  each  pod  wrapped  in  plastic  to  still  the  flow.  These  zero 
deployments  were  done  off  the  dock  on  March  19,  1989  after  the  experiment 
The  maximum  magnitudes  of  the  aero  velocities  are  3.77  cm/s.  One  channel 
from  each  pod  may  be  discarded.  A  large  offset  zero  is  an  indicator  of  trouble 
and  is  usually  the  channel  which  is  discarded.  Thus  die  largest  offset  zero 
which  is  used  is  -0.79  cm/s,  (pod  3  is  generally  suspect  on  both  tripods  and 
will  be  excluded  from  the  analysis  whenever  possible).  By  using  these  zeros, 
velocity  vectors  are  obtained  which  yield  logarithmic  profile  slope  with  regres¬ 
sion  coefficients  in  excess  of  0.98  indicating  confidence  of  better  than  0.2  cm/s 
in  the  mean  speed  estimates. 

Table  1 

Zeros  (cm/sec)  for  T2000  and  T2006  March  19  20:30  GMT 
pod  # 

channel  1  2  3  4  5  6 

A  1.4136  -0.1026  -0.0342  -0.6118  0.4370  -0.7638 

B  -0.4522  -0.1254  0.0760  -0.4028  0.7410  -0.1520 

C  0.7904  -1.2768  -2.6828  0.9500  0.6194  0.6270 

D  -0.8816  0.1406  -2.8614  0.4636  -0.7486  0.8816 

Zeros  for  T2005  (cm/sec)  March  19  18:30  GMT 
pod  # 

channel  1  2  3  4  5  6 

A  “  1.2388  3.7772  -0.8930  0.4712  " 

B  *  0.7714  -1.5352  0.7980  -1.2160  " 

C  -  0.4180  3.2414  -0.7828  “ 
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D  "  0.6270  0.2888  0.5130  0.1672  ' 

Figure  1  shows  a  few  sample  profiles.  The  four  symbols  for  each  level 
represent  speed  estimates  based  on  three  axes,  giving  four  combinations.  The 
quality  of  die  zero  calibrations  is  visually  represented  by  the  consistency 
between  estimates  of  speed. 


Transmissometers 

The  transmissometers  directly  measure  percent  transmission,  %Tr 
100%*bits/32768.  light  attenuation,  ctsm"1,  is  derived  from: 

m—OJ 


Tr 


Tr0e 


Tr . 


a  =  -“log”””* 
l  Tr 

where  /  is  the  path  length  of  tite  light  beam  and  Tr0  is  a  calibration  transmis¬ 
sion  at  zero  light  attenuation.  Tra  is  due  to  dear  water  attenuation  and  an  unk¬ 
nown  effect  from  plastic  spacers  we  used  to  reduce  the  pathlength  from  25  cm 
to  10  cm.  Their  attenuation  effect  has  not  yet  been  calibrated.  But  attenuation 
is  simply  offset  a  constant  amount  by  Tr0.  The  preliminary  plots  use 
7ro*100%  Therefore  when  die  transmissometers  were  working  the  measures 
of  light  attenuation  can  be  seen  to  go  from  a  minimum  of  "4  m-1  which  is  the 
offset  due  to  Tr0  to  maximum  values  of  30  m_1  which  represent  light  attenua¬ 
tion  of  at  least  26  m-1. 


Thermistors 

Thermistors  will  be  properly  calibrated  later  this  FalL  A  rough  calibration 
has  been  applied: 

T  =  25*Wts/65535 

where  bits  ranges  0  to  65535  corresponding  to  0  to  25  degrees  Centigrade. 


Pressure  FFT  calibrations 


The  Fourier  transform  of  the  pressure  signal  is  calibrated  with  conversions 
of  measured  bits  to  meters  of  salt  water  and  die  various  scaling  factors  con¬ 
tained  in  the  Fourier  transform  and  its  gain  adjustments.  Ultimately  die  sum¬ 
mation  over  die  components  of  a  Fourier  transform  give  die  variance  of  the  sig¬ 
nal: 


**  freq*Q 


Where  P2  are  the  squared  magnitudes  of  die  pressure  spectra.  These  are  writ¬ 
ten  to  tape  in  minimal  form  and  are  called  nff tbits2*.  The  gain  and  FFT  cali¬ 
brations  reduce  to: 


c}(bits)  =  -r—  £//fMtt2/20 
**fnq 

The  final  conversion  is  from  bits  to  meters  of  salt  water.  The  ParoScientific 
pressure  signal  is  a  count  given  as  a  3  byte  number,  C  The  count  is  converted 
into  a  period,  T,  by  inverting  the  frequency  found  by  the  product  of  count  and 
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sampling  rate: 

T*— rr“ — r  microseconds 

( C*rate ) 

The  pressure  is  given  by  the  equation: 

p  =A(i-TcrryB(i-T0rD2 

Where  A,  B  and  7*  are  calibration  coefficients  with  units  of  PSIA,  PSIA  and 
microseconds,  respectively.  The  conversion  of  PSIA  to  meters  of  salt  is 
achieved  with:  68.947  millibar  «  1.  psi;  1  millibar  *  (1.0197)"1  cm  water  or  1 
millibar  ■  (1.0197*  1.02813)"1  cm  standard  salt  water,  1  atmosphere  ■  10133 
millibar.  So: 

p  *  (P*68.947-10133y(1.02813*  1.0197)  cm 

A  few  pressure  numbers  in  counts  from  data  tape  t2006  where  die  depth 
should  be  approximately  90  meters  are  (01647f,  01646c,  01648b)bex  *  (91263, 
91244,  91275)decimaL  If  rate  ■  .4  these  numbers  give  about  73  rotters.  The 
slope  of  die  curve  is  -47.636  bits/meter  salt  water.  Hg  1  is  the  calibration 
curve  of  meters  salt  water  to  counts  recorded,  or  bits.  Uncertainly  in  this  cali¬ 
bration  is  due  to  my  most  recent  calibration  factors  for  the  pressure  sensors 
being  from  1982. 

Although  die  full  calibration  for  ParoSciendfic  pressure  sensors  is  non¬ 
linear  the  range  of  fluctuations  about  die  mean  of  "90  meters  is  very  nearly 
linear.  The  slope  is  very  near  -50  bits/meter  salt  water.  The  full  calibration 
conversion  is 

afQneters 2)  ■  £  iff  tbits  /50bia  /meter  )*/20 

ttfrtq 

Derived  quantities: 

The  four  components  of  die  velocity  vector  measured  by  a  BASS  pod  are 
rotated  into  three  component  Cartesian  velocity  vectors: 

*  A* 


where  A  is  die  rotation  matrix  which  takes  the  four  components  of  a  BASS  pod 
into  the  three  Cartesian  components.  A  may  be  defined  five  ways.  Four 
matrices  are  based  on  combinations  of  just  three  axes  at  a  time.  The  fifth  is  die 
average  of  die  ocher  four  methods: 

1  0  -1  o' 
a54  o  -i  o  i 

^  [03  03  03  03 

A  rotation  to  put  die  U,V,W  vector  into  a  NorthJsast,Up  coordinate  sys¬ 
tem  is  based  on  die  compass  reading.  The  compass  is  measured  in  the  pressure 
case  of  the  BASS.  The  velocity  vector  begins  in  coordinates  with  U  aligned 
with  die  A-C  axis.  The  angle  of  the  A-C  axis  with  the  compass  case  for  the 
Clean  tripod  is  180°  +  15°.  U  will  be  turned  into  East  so  that  is  another  90°. 
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Magnetic  declination  at  the  STRESS  area  is  -16.9*  in  1988-9.  The  compass 
angle  measured  for  T2006  was  266*.  The  total  rotation  is  therefore  8  ■  90®  - 
(266*  +195° )  - 16.9*  ■  332*.  This  angle  rotation  is  around  die  Z  axis: 
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The  ten  components  of  the  cross  correlations  of  the  four  axis  of  the  BASS 

pod  are  rotated  into  die  six  unique  components  of  the  Reynolds  stress  tensor. 
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Futber  rotations  are  applied  to  put  iqify  into  streamline  coordinates. 

A  least  squares  fit  of  die  logarithmic  velocity  profile  equation. 
U(z)  m  iutK  ln(z/z0),  yields  the  shear  stress  velocity  scale,  u*,  the  roughness 
length  scale,  z9%  and  the  regression  coefficient,  R2.  Only  the  current  meters 
below  2  meters  were  used  as  the  velocity  profile  deviates  from  logarithmic 
above  that  level. 


Linear  wave  theory  can  be  applied  to  solve  for  the  near  bed  wave  orbital 
velocity.  ub,  near  bed  wave  induced  panicle  excursion  distance,  ab,  and  the 
surface  wave  amplitude,  h,  given  die  near  bottom  pressure  variance,  Oj>  (meters 
«  P/  p  g),  and  the  wave  frequency,  f,  obtained  from  the  spectra. 

k  * 


<o  =  2 a/ 

.  2o.  « 

b  r  sinhikh) 

ah  k  2o»  3&VSL 

b  f  sinh(kb) 
h  =  2jof  cosh  (kb) 


The  wave  velocity  ub  will  be  the  major  contribution  to  the  kinetic  energy  or 
die  trace  of  the  Reynolds  stress  tensor.  If  this  is  true  then 
,2  _  —  J_ 


q2  =  uV+vV+wV  =  -fa*2. 
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Fig.  1.  Sample  speed  profiles  showing  the  four  possible  calculations  of  speed  based  on  choosing 
three  of  the  four  BASS  axes,  (symbols  o,*,x,+).  When  one  axis  is  bad,  three 
methods  of  calculation  use  that  axis  and  result  in  erroneous  estimates. 
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Fig.  6.  Two  minute  average  velocity  at  1.1  mab  (solid  lines),  thirty  minute  averages  (o),  forty 
hour  low  pass  filtered  (+).  This  sample  shows  dramatic  internal  waves  when  the 
mean  velocity  was  small. 
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logarithmic  velocity  profile  yield  u.  and  z„.  This  plot  shows 
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11  Abend  (IMC  S00  worde) 

Two  Benthic  Acoustic  Stress  Sensor  (BASS)  equipped  tripods  were  deployed  in  the  Sediment  Transport  Events  on  Shelves 
and  Slopes  (STRESS)  experiment  in  November,  1988,  and  recovered  in  March,  1989,  on  the  California  Shelf.  They  measured 
velocity  profiles  in  the  bottom  boundary  layer  over  the  lowest  3  meters.  Transmissometers,  thermistors,  and  a  pressure  sensor  on 
each  tnpod  provided  suspended  sediment  concentration,  stratification,  and  wave  spectral  information,  as  well. 
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